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Abstract 

We construct solvable models on the honeycomb lattice by com- 
bining three faces of the square lattice solvable models into a hexagon 
face. These models contain two independent, anisotropy controlling, 
spectral parameters and their transfer matrices with different spectral 
parameters commute with each other. At the critical point, the finite- 
size scaling of the transfer matrix spectrum for the honeycomb lattice 
models is written in terms of the quantities obtained from the finite- 
size scaling of the square lattice solvable models. We study in detail 
the phase transition properties of two models based on the interacting 
hard square model and the magnetic hard square model, respectively. 
The models, in general, can be extended to the IRF version of the 
Z-invariant models of Baxter. 



I Introduction 



The Yang-Baxter equation, also called the star-triangle relation (STR) has 
been one of the most useful tools for analytical calculation in the statistical 
lattice models|l|. A lattice model is called solvable if its Boltzmann weights 
satisfy the STR. There are a variety of solvable models. These models are all 
defined on the square lattice. In this work, we construct solvable models on 
the honeycomb lattice by combining three faces of the square lattice solvable 
models into a hexagon. The working is best illustrated from the perspective 
of the kagome lattice eight- vertex model and its Ising spin version on its dual 
lattice as follows. 

In 1978, Baxter suggested the Z-invariant model on an arbitrary planar 
lattice where any three lines do not intersect at a point. || The intersections 
of the lines are regarded as lattice sites and an arrow is assigned to each 
section of line divided by two other lines. A model is said to be Z-invariant 
if its interaction satisfies the STR of vertex model for any triangle of the 
planar lattice. Here Z-invariant implies that the partition function Z remains 
unchanged after sliding a line over a vertex point. One of the examples of Z- 
invariant model is the kagome lattice eight- vertex model which is investigated 
in Ref.0. 

On the kagome lattice, there exist three types of vertices. For solvable 
eight-vertex model on the kagome lattice, each type of vertices is described 
by different spectral parameters and is related to other types by the STR. 
The total free energy of this model is simple sum of three free energies each 
of which are obtained from the regular square lattice eight- vertex model of 
the corresponding spectral parameter. The two-point correlation function 
(cr r a m ) of two spins a r and a m located at site r and m depends only on one 
spectral parameter, if the sites r and m lie between two adjacent parrallel 
lines and the system is large enough. 

The kagome lattice eight-vertex model can be transformed to the inter- 
action-round-a-face (IRF) model on its dual lattice. The dual lattice sites 
can be divided into two kinds of sites; one on the face of the hexagons and 
the up triangles and the other on the face of the down triangles of the kagome 
lattice. The former forms the honeycomb sublattice and the latter the tri- 
angular sublattice. The dual lattice of the kagome lattice will be called the 
honeycomb-triangular lattice. We can assign Ising spins to the honeycomb- 
triangular lattice points according to the direction of arrows of the original 
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vertex model. Since the vertex weight of the eight-vertex model is deter- 
mined by four arrows surrounding a lattice point, the corresponding Ising 
model possesses both two and four spin interactions. When one takes the 
partial trace over the spins on the triangular sublattice, the resulting system 
is an Ising model on the remaining honeycomb lattice with complicated in- 
teractions. In this way, one can obtain the honeycomb lattice Ising model in 
which a hexagon is the unit of the IRF model. 

We write the face weight of the Ising model corresponding to a vertex 
with spectral parameter u as w(a, b, c, d\u) where a, b, c and d denote the spin 
states located at lower-left, lower-right, upper- right, and upper-left corners 
of a face, respectively, and denote by W^a, b, c; d, e, /) the hexagonal face 
weight for a face whose corners have the spin values, a, b, c, d, e and /, 
counterclockwisely. Then the face weight for a hexagonal face is written in 
terms of the three spectral parameters u, v and z as (See Section II and 
Fig.0.) 

W h (a, b, c; d, e, /) = ^ w(e, g, c, d\u)w(g, a, 6, c\v)w(f, a, g, e\z) . (1) 

9 

The summation runs over the possible spin values ±1. Following the condi- 
tion for the Z-invariance JTJ , the spectral parameter z is connected to the other 
parameters u and v by z = u — v . In the notation of W^, two independent 
variables u and v will be specified as W^a, b, c; d, e, f\u, v). 

An important observation is that the model obtained in this way is also 
solvable as we discuss in Section II. In fact, one can construct honeycomb 
lattice solvable models (HLSM) in this way starting from arbitrary weights 
w(a, b, c, d\u) of solvable square lattice models. The purpose of this paper is 
to formulate generally the HLSM and study their properties. This paper is 
organized as follows. In Section II, we define the HLSM in a general setting 
and show the commutivity of the row-to-row transfer matrix. We also derive 
the general form for the bulk and the finite-size corrections of the transfer 
matrix spectra. In Section III and IV, we study the HLSM's based on the 
interacting hard square and the magnetic hard square, respectively. The 
ground state symmetries are mapped in the (u-v) plane and, in the case of 
the magnetic hard squares, the transfer matrix spectrum is analyzed. In 
Section V, we extend the HLSM to a hierarchy of models which have larger 
unit cells. We conclude with a discussion in Section VI. 
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Figure 1: Graphical representation of the star-triangle relation. 

II Honeycomb lattice solvable models 

We start with the definition of the honeycomb lattice solvable models (HLSM) 
which are the IRF models constructed by using solvable models on the square 
lattice. Let (a, b, c, d) denotes a spin configuration of a unit square for a lat- 
tice spin model on the square lattice. Values of a, b, c, d range over the set 
of spin states depending on models. The face weights of solvable IRF models 
are denoted as to (a, b, c,d\u) where u is the spectral parameter and its depen- 
dences on other interaction parameters are implicit, w is the solution of the 
STR which is shown graphically in Figfj], and satisfy the standard reflection 
symmetries 

w(a, b, c, d\u) = w(c, b, a, d\u) = w(a, d, c, b\u) 
and the initial condition 

w(a,b,c,d\0) = 5 ac . 

For any solution w(a, b, c, d\u) obtained from the STR, we can construct the 
HLSM whose hexagonal face weights are given by Eq.fjj]). This definition 
of the hexagonal face weights is shown diagramatically in Fig.[| Since the 
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Figure 2: Definition of the Boltzmann face weights in the honeycomb lattice 
solvable model. 

square lattice face weights satisfy the STR, we get the identity of 

W h (a, b, c; d, e, f\u, v ) 

= a ' & ' 9\u)w(g, 6, c, d\u - v)w(e, /, g, d\v) . ( 2 ) 

9 

Using the diagonal reflection symmetry to Eq.(|2|), we have the symmetry 
relation 

Wh(a, b, c; d, e, f\u, v) = Wh(d, e, /; a, b, c\u, v) . 

We consider a honeycomb lattice system with M rows each of which is 
made of iV hexagons as shown in Fig.|| To the lower and the upper sites of 
a row, we assign the spin configurations a = (ai, CL2, • ■ • , a2v) and a' = (a[, 
a' 2 , • ■ • , a' 2N ), respectively. The a^'s and a^'s are the possible spin states. For 
convenience, we count the sites so that the set of spins {a 3 , 04, as, a' 4 , a' 3 , 
a' 2 } constitutes a hexagonal face. In this paper, we restrict ourselves to the 
periodic boundary condition; a 2 7v + i = ai and a' 2N+1 = a[. 

Defining the transfer matrix V^(«, v) as 

V SS /(w,u)= Yl w h(ai,ai+i,ai+2]a' i+1 ,a' i ,a' i _ 1 \u,v) , 

Odd i<2N 

the partition function Z of the system with M rows is written in terms of 
the transfer matrix V(w, v) as 



Z = tr[V(u,v)] 



M 
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Figure 3: Distribution of the spins to the honeycomb lattice sites. a x: a t 
and #h are the parameters appearing in the finite-size scaling and not to be 
confused with the spin state labelings. 

The transfer matrix V of the honeycomb lattice model is the function of 
two independent spectral variables u and v. Using the STR, the commuting 
relation of 

[V(u,v),V(u',v)] = 

can be proved for any values of u and it'. This is best illustrated graphically 
as shown in Fig.§. In a similar manner, we also get the commuting relation 
of 

[V(u,v),V(u,v')] = 

for different v and v '. Those relations imply the integrability of the model. 
The free energy per a hexagon face is given by 

/h(«, v) = f(u) + f(v) + f(u - v) , 

where /(it) is the free energy per square of the square lattice system defined 
by the original face weight w(a,b,c,d\u). This relation is derived in the 
same way as how the free energy of the kagome lattice eight-vertex model is 
obtained in [I]. 

Assuming the honeycomb lattice model has the conformal symmetry at 
the critical point, the finite-size scaling of the transfer matrix spectrum can 
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Figure 4: Proof of the relation \V(u, v), *V(u', v)] = 0. (a) Two rows with 
the transfer matrix *V(u,v) and V(u',v), respectively, (b) By the STR, we 
have moved the faces with the spectral parameter v over the faces with the 
parameter u and u — v . (c) Two rows in the middle have been interchanged, 
using the properties of square lattice solvable models, (d) By the STR, the 
uppermost faces with the spectral parameter v go down beyond the faces 
with the spectral parameter u' and u' — v. 
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be predicted by the conformal field theory. For the a-th eigenvalue A„ of the 
transfer matrix V(u,v), we expect the finite-size scaling form as follows 0, 

-\n\A a \=Nf h + ^(x a -^)^sm9 h (3) 

under the periodic boundary condition. Here /h is the bulk free energy, c 
is the central charge and x a is the scaling dimension associated with the 
corresponding correlation function. a x is the width of one hexagon in the 
horizontal direction and at is the length in the time direction. a x and at 
are illustrated in Fig|| 9^ is the angle depending on the anisotropy of the 
interaction 0]. For general u and v, the lattice should be deformed to make 
the system isotropic in the continuum limit so that a t , a x and 9^ take non- 
trivial values. 

In Ref. 0, it was shown that each unit square of solvable square lattice 
models is deformed in such a way that a t = a x due to the diagonal reflec- 
tion symmetry and the deformation is characterized by the anisotropy angle 
9 which depends linearly on the spectral parameter. Applying this result 
to each types of squares forming a hexagon, one finds that the anisotropy 
parameters in Eq.@ for the HLSM are given by 

a T = la cos ■ 



a t = 2a cos 
9 6{U) 



2 

9{u-v) 



where 9(u) is the anisotropy angle of the square lattice solvable model with 
the spectral parameter u and a is the lattice constant of a hexagon, that 
is, the distance between nearest neighbors of the honeycomb lattice. These 
relations yield the finite-size scaling of the honeycomb lattice model, in terms 
of the quantities obtained in the square lattice solvable model as 

-ln|A«| = N{f(u) + f(v) + f(u~v)} 

2tt , c .cos{9(u — v) /2} . . , 

+ N - 12» L Wl ,)/2 } Sm ^' V > ■ < 4 » 

Due to the initial condition w(a, b, c, d|0) = 5 ac , the limit of v — > reduces 
the honeycomb lattice to the square lattice system of the spectral parameter 
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u with 2N faces in a row. In the limit v = 0, the above finite-size scaling 
relation Eq.(||) reduces to that of the square lattice model with the parameter 
u as it should be. 

One of the simplest examples for the honeycomb lattice solvable models is 
the standard Ising model. The solvable Ising model with diagonal interaction 
on the square lattice is defined by the interactions K u and L u through the 
face weight 

w(a, b, c, d\u) = exp(K u ac + L u bd) , 

where the ^-dependences of K u and L u are given in the standard references, 
e.g. Ref.0. For the independent parameters u and v, the hexagonal weight 
is defined as 

W h (a,b,c; d,e,f\u,v) 
= exp(K u ec + K u „ v fg + K v gb + L u gd + L u _ v ea + L v ac) . 

9 

The STR simplifies this weight to 

Wh(a, b, c; d, e, f\u, v) = exp(K u ce + L u _ v ea + L v ac) 

x exp(K u fb + L u _ v bd + L v df) 

which shows that only the two body interactions between next nearest neigh- 
bors appear on the honeycomb lattice. Since there is no interaction between 
nearest neighbors, the honeycomb lattice Ising model is equivalent to two in- 
dependent Ising models on the triangular lattice. However, for other solvable 
models, such simple decomposition does not arise and following two sections 
are devoted to specific models, i.e., the interacting hard square model and 
the magnetic hard square model. 

Ill Interacting hard square (IHS) model 

In this section, we consider the HLSM based on the interacting hard square 
(IHS) model which is a diagonally interacting lattice gas model where a site 
cannot be occupied unless its nearest neighbors are all empty. Each site has 
the value O(empty) or 1 (occupied) and the face weight is described by the 
fugacity z, the interaction L and M: 

w(a, b, c, d) 
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f mz (a+b+c+d)/4 exp (£ ac + Mbd)t- a+b ~ c+d , if ab = be = cd = da = 0, 
| , otherwise. 

Here m is a constant which can be adjusted for a suitable normalization. 
The quantity t is necessary so as to discover the solution satisfying the STR. 
In the homogeneous square lattice where all the square faces have the same 
interaction strength, t does not appear in the partition function. However, 
on the honeycomb lattice, it can play a non-trivial role. Baxter[|J found the 
solutions to the STR on a special manifold in the parameter space (z, L, M, t). 
The manifold is parametrized by the temperature-like variable q and the 
spectral parameter u. Interaction parameters, hence, the face weights are 
expressed by elliptic functions with argument u and nome q 2 . The face 
weights are defined when the spectral parameter, u varies in the range of 
— 7r/5 < u < 27r/5. 

On the honeycomb lattice we will label six sites of a unit cell, counter- 
clockwisely from the uppermost corner, by sequential numbers 1 to 6. For 
instance, in Fig.|| we assign the site of the spin state e as 1, / as 2, and so 
on. And we denote by the interaction between the sites i and j. is 
the three-spin interaction among the sites i, j and k. Then the face weights 
for the honeycomb lattice IHS model are given as 

W h (a,b, c; d,e,f\u,v) 

(a+b+c+d+e+f)/3 i T 1 , t £ i r J\ 

= mo^h exp(L 14: eb + L 2 5jc + L 36 ad) 

x exp(Li 3 ea + L 2i fb + L 35 ac + L m bd + L 51 ce + L 62 df) 
x exp(Li 35 ace + L 2m bdf) 

where Zh is the fugacity and mo is the weight for the vacuum state, i.e. the 
state where all six sites of a hexagon are vacant. Actually the vacuum state 
of this honeycomb lattice model contains the state with one particle at the 
center of the cell. For the solvable IHS model, there is no four-spin interaction 
on the honeycomb lattice. From the exact solution of the square lattice, we 
obtain the following functional form for the interactions: 

- u)d 1 {u - v)d 1 {v) 
+ 0?(f)0i(¥ +u)*i(£ + u-v)9 1 Cf + v) 
O^WuWi* - u + v)%(z - v) 3 
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-m 



where Q\ is the standard elliptic theta function. It is defined by the spectral 
parameter u and the temperature-like parameter q 2 as 

00 

0i (u) = sinw -V n cos2u + g 4n )(l -g 2n ) . 

ra=l 

For simplicity, we shall omit q 2 in the notation of Q\. It is to be noted that 
the interactions show the symmetry under the rotation by angle 180° which 
leads to the property of Wh{a, b, c; d, e, f\u, v) = Wh(d, e, /; a, b, c\u, v). 

q 2 — > — 1 is the high temperature limit and q 2 = 1 corresponds to zero 
temperature. If we compute the hexagonal face weights for any u and v in 
these two limits, the dominant configurations on the honeycomb lattice can 
be obtained. By tracing the dominant configurations, we can draw the phase 
diagram of honeycomb lattice model. There are four regions in the (u-v) 
plane where the ground state configurations take different nature. This is 
shown in Fig.|5|. The ground state configurations of each region are drawn 
in Fig|| The region D includes the point of (u,v) = (2n/15, 7r/15) where 
the interaction becomes isotropic. Fig.|5](D) shows a ground state of the 



11 




Figure 5: Phase diagram of the honeycomb hard square model. A, B, C and 
D are the regions in which the exact solution of the IHS model is defined. 
In order to emphasize the three-fold symmetry of the lattice, we choose the 
non- rectangular coordinate system where f-axis makes an angle 120° with 
the w-axis. 
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Figure 6: The ground states in the interacting hard square model. The figure 
(A), (B), (C) and (D) are the ground states of region A, B, C and D, of Fig.5, 
respectively. 
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region D. Even if we interchange the empty and occupied sites in the state 
of Fig.|||(D), it does not affect on the energy of the system. That is, there are 
doubly degenerate ground states. Their particle density is 1/2. The critical 
point, corresponding to q = 0, in this region is of the tricritical Ising model 
universality class as in the original square lattice model. In the regions A, 
B, and C, the particle density of the ground state is 1/3. The ground states 
of the region B and C are obtained by rotating the system of region A by 
120° and —120°, respectively. The honeycomb lattice is divided by three 
sublattices and the ground states of the region A, B, and C have three-fold 
degeneracy. As the temperature increases from zero, the domain walls will 
appear. The system having three different domains is described by the 3- 
state Potts model. In the high temperature limit of q 2 — > — 1, we can see 
as the dominant phase the vacuum phase in the regions A, B and C. 

Let us call the ground states of (A), (B) and (C) as 3 x 1 ordering and 
that of (D) as v2 x a/2 ordering. We borrow this notation from the square 
lattice ordering. When v approaches to 0, the system becomes the square 
lattice and the ground states given in Fig.| are just the 3x1 and y/2 x y/2 
orderings on the square lattice. 

The 3x1 orderings on the honeycomb lattice may look rather strange. 
However, by converting the anisotropy of interaction to the geometrical 
anisotropy of lattice as discussed in Ref.[[J, they are reduced to the com- 
pact ordering of large particles as shown in Fig.[7[ Therefore, the nature of 
symmetry breaking in the honeycomb lattice is the same as in the square 
lattice and no new universality classes appear in the HLSM. 



IV Magnetic hard square (MHS) model 

The magnetic hard square (MHS) model is the IHS model where each of 
particle has spin 1/2 internal degrees of freedom. When a site is empty, its 
spin value is assigned to be zero. Otherwise, the spin is 1 or —1. Its face 
weights are defined by 

w(a, b, c, d) 

f mz (* 2 +b 2 +c 2 +d*)/4 eX p( La 2 c 2 + Mb 2 d 2 + Jac + Kbd)t^- b2+c2 - d ^l A 

— I , if ab = be = cd = da = 0, 

[ , otherwise. 
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Figure 8: Definition of four regions A, B, C and D in the honeycomb lattice 
MHS model 

The exact solution of this model on the square lattice which obeys the STR 
was found by Pearce||. 

For the solution on the so-called E manifold, the face weights are written 
by the elliptic theta functions parametrized by the spectral parameter u and 
the temperature-like parameter q. On the E manifold, the spectral parameter 
varies in —tt/6 < u < tt/2. q has the value between —1 and 1. On the 
square lattice, q — > 1 is the zero temperature limit and q — > — 1 is the high 
temperature limit. At q = 0, the E manifold crosses the critical line of T 
manifold where all the face weights are the trigonometric functions of u. The 
T manifold is a line parametrized by the crossing parameter A and the q = 
point of the E manifold corresponds to the A = 7r/3 point of the T manifold. 

Applying the elliptic solution of the square lattice to the honeycomb lat- 
tice model, we can draw the phase diagram of the honeycomb lattice MHS 
model. It is divided into four regions as shown in Fig|| At zero tempera- 
ture, the particles form \/2 x \/2 ordering as in the previous section. In the 
\/2 x a/2 ordering the honeycomb lattice has two sublattices. The ground 
states are that the sites of one sublattice are occupied and those of the other 
are empty. In the region D including the isotropic point, the system shows 
the ferromagnetic phase at zero temperature limit. In the other regions the 
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(A) (B) (C) 



Figure 9: Antiferro magnetic orderings which are the ground states of the 
region A, B and C, respectively. 

antiferromagnetic phases are dominant where +1 and —1 spins are alternat- 
ing line by line. In the square lattice limit, this corresponds to the (0101) 



state, following the notation used in Ref.JTT]]. Fig.|9| illustrates the antiferro- 
magnetic phases of the MHS model. In the high temperature limit, the region 
D is the paramagnetic phase where the spins of the occupied sublattice are 
distributed randomly and the local magnetization is zero. This limit belongs 
to the non-magnetic IHS model discussed in the previous section. Since there 
is no magnetic interaction in the IHS model, neither spin is preferred and 
the system shows the random distribution of the spin. In the other regions, 
the vacuum state is dominant. 

The T manifold is the critical line where the face weights are the trigono- 
metric functions of the spectral parameter. On the T manifold, the bulk free 
energy is given by that of the six-vertex model. This manifold is classified 
according to the spectral parameter u into three regions called as I /II, III/IV 
and V/VI region.^ In I/II region, u varies in — (n — A) < u < and in III/IV 
region, in < u < A. V/VI region is where u is in the range A < u < tc. In 
each region, the anisotropy angle of the square lattice which appears in the 



lr The word regime is conventionally used instead region for classifying the regions ac- 
cording to the spectral parameter. But the word region would be more proper to imply 
the range depending on the spectral parameter. In this paper, we choose the word region. 
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finite-size scaling is found [[| ||, [| to be 








in region I/II, 
in region III/IV, 
in region V/VI. 



(5) 



The central charge of the square lattice model is c 



10- 



By calculating numerically the eigenvalues of the transfer matrix for the 
honeycomb lattice model and using the finite-size scaling proposed in Sec. II, 
the central charge and the scaling dimensions are obtained. From the nu- 
merical results, this model is manifested to have the central charge c = 1. 
In the case of A = 5n/12, u = A/3 and v = A/8, we calculated the eigenval- 
ues for N = 4, 5 and 6. It is the case where all of three faces constituting 
one hexagonal face are in III/IV region. Table [I] shows the data of iV = 6 
for the scaling dimensions and the extrapolated values. We extrapolate the 
results by fitting finite size data to the form x a (N) = x a (oo) + a/N b with 
a constant a and b to get the N- independent scaling dimension x a (oo). We 
divide the levels into sectors denoted as [S,f]. S is the spin of the field and 
f is the eigenvalue for the spin reversal operator. Their rigorous definition 
is given in Ref. 0. The starred levels in the table represent doublets. For 
other notations, see Ref . and |§ . The transfer matrix spectrum shows that 
the honeycomb lattice MHS model belongs to Ashkin- Teller class as in the 
case of the square lattice model. The last column of the table provides the 
corresponding scaling dimensions of the Ashkin- Teller model. (0, 0) is the 
level related to the identity operator. X„ jTO 's are the scaling dimensions of 
the gaussian operators, which are given by 



where n, m are integers and g is the gaussian model coupling constant and is 
related to the crossing parameter by g = 9(1 — A/71"). The dimensions (1/16, 
1/16), (9/16, 9/16) and (25/16, 9/16) are associated with the magnetization 
operator and are A-independent. 

When all the values of u, v and (u — v) are in the region I/II or V/VI, 
we also have the same result as above. But in the cases that the faces of 
the region III/IV are mixed with those of the region I/II or V/VI, it appears 
that there exist an additional factor in the bulk free energy term besides 
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[S,r] 


N = 6 


Extrapolated 


Exact value 


[0,+] 


-1.6 x 1(T 4 


1.4 x 10~ 5 





(0,0) 




0.0953 


0.0952 


0.0952 






0.3803 


0.3810 


0.3810 


X 2 ,0 




0.8725 


0.8543 


0.857 


-^3,0 




1.501 


1.523 


1.524 


-^4,0 




1.965 


2.000 


2 


(1,1) 




2.040 


2.099 


2.095 


*i,o + 1 + 1 




2.317 


2.374 


2.381 


-^5,0 




2.342 


2.373 


2.381 


X 2 ,o + 1 + 1 




2.562 


2.620 


2.857 


X 3 , + 1 + 1 


[o,-] 


*0.1254 


0.1251 


0.125 


(1/16,1/16) 




*1.127 


1.128 


1.125 


(9/16,9/16) 




♦2.121 


2.131 


2.125 


(1/16 + 1,1/16 + 1) 


[1,+] 


1.089 


1.094 


1.095 


x lfl + 1 




1.362 


1.377 


1.381 


X 2 ,o + 1 




1.861 


1.870 


1.857 


X 3 , + 1 




2.421 


2.503 


2.524 


X 4i0 + 1 




2.652 


2.706 






2.887 


2.955 






[1,-] 


*1.120 


1.125 


1.125 


(1/16 + 1,1/16) 




*2.069 


2.099 


2.125 


(9/16 + 1,9/16) 




*2.092 


2.131 


2.125 


(25/16,9/16) 


[2,+] 


1.962 




2 


(0 + 2,0) 




2.010 


2.078 


2.095 


X lfi + 2 




2.065 




2.095 


X lfi + 2 




2.222 


2.624 


2.381 


X 2 fi + 2 




2.343 




2.381 


X 2 ,o + 2 


[2,H 


*2.018 


2.112 


2.125 


(1/16 + 2,1/16) 




*2.080 


2.112 


2.125 


(1/16 + 2,1/16) 



Table 1: Transfer matrix spectrum of the honeycomb lattice MHS model 
with A = 5n/12, u = A/3 and v = A/8. The extrapolation is done from the 
data of iV = 4, 5, and 6. The starred levels are doublets. 
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the sum of the bulk free energy of each face. Since we could not calculate 
the spectra for large enough system, we could not resolve the nature of this 
anomaly. Resumably, it comes from the mismatch between the ground states 
of different regions. Even if there exist an unexplained additional term, the 
1/N correction term in the finite-size scaling show the expected behavior of 
Ashkin- Teller class for all the cases we have calculated. 

The MHS model is generalized to the affine D model. The affine D model 
is characterized by A and an integer L. The number of spin states is (L + 3). 
The MHS model is equivalent to the affine D model with L = 3. In the general 
affine D model, the exact solution is found by Kuniba and Yajima|TT|] when 
A = Ti /L. (Pasquier's ADE model corresponds to A = 0.[|K|) The transfer 
matrix spectra of the square lattice affine D model are analyzed in Ref. [|T^, §]. 
It is found that they belong to the Ashkin- Teller universality class with the 
gaussian coupling given by g = L 2 (l — X/tt). For the general affine D model, 
the honeycomb lattice model is also defined. We have also analyzed the 
transfer matrix spectrum for several cases of the honeycomb lattice affine D 
model and have found that their critical properties are again the same as in 
the square lattice model. 



V Extension of the honeycomb lattice solv- 
able models 

Up to now we have investigated some properties of the honeycomb lattice 
models where each hexagonal face weights are given by the multiplication of 
the square face weights with the parameter u, (u — v) and v. In this section 
we extend this scheme of constructing the integrable models to another kind 
of integrable models with larger unit cells. A hexagonal face consists of three 
parallelograms and we will distinguish the spectral parameters for the upper- 
right, upper-left and lower parallelograms by u, u' and u , respectively. For 
the models defined in Sec. II, they are fixed to be u' = u — v and u = v. 
When the size of a unit cell is given as m x n with positive integers m and 
n, let us label the spectral parameters of the hexagon located at z-th column 
and j-th row of a unit cell as Uij, u'^ and u^. Here is the positive integer 



modulo m(n). In Fig.|10|, the unit cell of 3 x 4 hexagons is illustrated as an 



example. Suppose that the spectral parameters of parallelograms included 
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Figure 10: A unit cell consisting of 3 x 4 hexagons 

in the unit cell are related by two variables of (u, v) as following: 

Uij = u + (i — 1) • (u — v) 

u'ij = -(n-j + l)v (6) 
u"j = (n-j + l)v . 

As long as the spectral parameters satisfy Eq.(|6]), the integrability of this 
model also can be proved by the STR in the simple way Thus one can 
generate a hierarchy of integrable models starting from each integrable square 
lattice model. 

As the honeycomb lattice solvable model corresponds to the deformed 
honeycomb lattice as shown in Fig.[7|, the system whose unit cell is made 
ofmxn hexagons with different anisotropy angles results in the periodic 
deformation of honeycomb lattice. One example is illustrated in Fig.|lT] for 
the case m = 1 and n = 4. 
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Figure 11: An example of the integrable model with 1x4 unit cell where 
the anisotropic interaction is interpreted as the geometrical anisotropy. The 
thick line denotes boundaries of unit cells. 



22 



VI Discussion 



In this paper, we generalized the honeycomb lattice Ising model which is 
the IRF version of the kagome lattice eight-vertex model and have obtained 
the general honeycomb lattice solvable models. These models includes two 
spectral parameters which control deformations of a square into a parallel- 
ogram of general shape. The finite-size scaling behavior of the free energy 
and the anisotropy factor are written in terms of the quantities obtained in 
the square lattice models. From the several solvable models satisfying the 
STR, we analyzed the corresponding honeycomb lattice models. As to the 
critical phenomena, we do not expect new universality classes appear in the 
models constructed in this way, because the two point correlation function of 
the honeycomb lattice model behaves as in the regular square lattice model. 
However, some noble behaviors may appear in the non-universal quantities 
if the system has different properties in the different regions of the spectral 
parameter. 

The IRF version of the Z-invariant model is described by the planar lat- 
tice system where the parallelograms of arbitrary angles fill the plane with 
no overlap and no gap. Considering a hexagon in such a lattice system, we 
can change the positions of three parallelograms constituting the hexagon 
by applying the STR. This process does not change the partition function. 
It allows us to construct lots of different systems having the same partition 
function, that is, the same free energy. If we consider the IRF version of the 
kagome lattice model, there are three kinds of the parallelograms. Keeping 
the free energy constant, we can change the distribution of the parallelo- 
grams by using the STR, repeatedly. Configurations obtained in this way 



will be the configurations of the random tiling model discussed in Ref. 
The Boltzmann weights of solvable models are analogous to the chemical 
potentials of the tiles. It also can be interpreted as the state of the zero- 



temperature triangular Ising antiferromagnet flI3| , if we use only three types 



of parallelograms in the solvable models. Possible application of the STR to 
the random tiling problems is left for future study. 
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